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Abstract 

In the standard hot big bang model, the expansion of the early universe is given by the Friedmann 
equation with an energy density dominated by relativistic particles. Since in a variety of models 
this equation is altered, we introduce modifications in the Friedmann equation and show that we 
can constrain them using big bang nucleosynthesis data. When there is no neutrino/antineutrino 
asymmetry these modifications are tightly bounded but in presence of an asymmetry the bounds 
become much looser. As an illustration, we apply our results to a model where the second and 
third families couple to gravity differently than the first family (non- universal gravity). 
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I. INTRODUCTION AND CONCLUSIONS 



ds^ = de - a" it) \ ^ + de^ + r2 sin^ ed<P^\ (1) 



The framework of our present understanding of the evolution of the universe is based 
on the hot big bang model. Assuming the properties of homogeneity and isotropy of the 

universe, the space-time metric has to be the Robertson- Walker metric 

dr"^ 
1 — kr"^ 

Here, a{t) is the scale factor, and k = +1, or —1 corresponds to a closed, flat or open 
universe. The other fundamental ingredient for the dynamics of the universe are the Einstein 
equations. These equations, along with the metric (Q), and in a universe containing a perfect 
fluid of total energy density p and pressure p, imply the Friedmann equation (FE), 

SttG k 



H'-^-^^'i. (2) 

with H = a/a the Hubble expansion rate. One also deduces the energy- momentum conser- 
vation equation 

P = -?>{p + p)H (3) 

To have a predictive model we need, in addition to Eqs. Q and (jHl), to use other inputs, 
like the density of baryons, the density of dark matter and energy, etc. However, Eqs. Q 
and are very basic equations directly determining the evolution of the universe. In this 
sense, it is worth to test them as much as we can, since they could be altered by the presence 
of new physics ll, S S 0, S]- In principle, we do not expect a radical change from the form 
((21), but the possibility of some alteration is open. 

In this article we will concentrate on the FE specifically on whether the evolution 
of the universe at its earliest times is given exactly by the FE. We will constrain potential 
modifications to the FE in the period of big bang nucleosynthesis (BBN), since this is the 
earliest epoch where we have observational information. We will show that we can test this 
equation with the present available data. 

The period of BBN starts at T ~ 1 MeV with the decoupling of the weak interaction 
processes that allow n < — > p inter-conversion and ends with the nuclear fusion production 
of the light elements, at about T ~ 50 keV. 

It is an easy exercise to particularize the FE (0) in the radiation era, where BBN happens. 
For T <^ m^, the relevant energy density appearing in Q is given by 

P = Po = P^ + Pu.+ Pu^ + Pur + Pe (4) 



In (0} we have the contribution from photons, 

f <Pk E _ , vr^ 

" y (2^ exp(E/T) - 1 " ^ 30 

with the degrees of freedom g-y = 2 and = |/c|. We also have the contribution of the three 
neutrino species (neutrinos plus antineutrinos) and of electrons plus positrons 

_ r d^k E r d^k E 

J j2^^exp[{E-i,,)/T] + l^^' J J2^^ exp[{E + fi,)/T] + l 

with the degrees of freedom g^ = 1 and ge = 2, E = ^A;^ + mf, and we allow for a particle- 
antiparticle asymmetry by introducing a chemical potential /ij. 
We write the FE in the radiation dominated epoch as 

H = po in 

with Po given in Q and where we have neglected the fc-curvature term in (0). 

The article is organized as follows. First, in section |nj we add a new contribution to po 
and constrain it when the neutrino/antineutrino asymmetry vanishes, p^ = 0. We show that 
the new contribution is tightly constrained. In particular we follow the suggestion of the 
model presented in j^, and introduce modifications in G. More specifically, we let the second 
and third generation neutrinos to enter in Eq. ((7j) with a coupling G2 and G3, respectively, 
different from G. As a consequence of our general study, we can place strong constraints on 
this model. In section UTTl we allow for p^, and show that one may have larger departures 
rom Eq. (jZj). However, even in this general case, we will be able to rule out the model of 



a 



II. TESTING THE FRIEDMANN EQUATION WITH NONDEGENERATE BBN 

Modifications to the FE arise in a variety of contexts. For example, in brane- world models 
with extra dimensions one obviously finds departures from standard general relativity, and 
hence the FE is altered. One may have radical changes in this equation when compared to 
the standard equation ^ [l| but then it is difficult to reconcile the theoretical predictions 
with observations. Thus, the hard work is in the struggle to find scenarios where one gets 
an evolution equation that at first approximation is close to the FE, with modifications that 



can be tuned to be small 



How small they must be ? We will give an answer using certain 



type of alterations of the FE (|7j), and constraining them in the radiation period using BBN. 
In this section we set fi^i = 0, but in the next section we will relax this condition. 

Apart from brane-world models, we have also found inspiration in two other recent papers, 
where the general FE (j21) is assumed to be valid, but where there are contributions to p not 
having exactly the form (0), valid for radiation. First, in ^, the authors show that the 
reheating temperature Tr can be as low as Tr ~ 0.7 MeV (when BBN starts). For T above 
Tr we would have the decay of coherent oscillations of the inflaton, or the decay of a modulus 
or similar field. If indeed happens to be as low as speculated in jj] , it is conceivable that in 
the period of BBN, for T < T^, apart from having radiation dominating the energy density p, 
we might still have some traces of the period T > amounting to some small contribution 
to p that should be added to the r.h.s. of Eq. ((Tj). Second, in another context, the so-called 
Cardassian expansion model has been shown to lead to a model in which the universe is 
fiat, matter dominated and accelerating Qj. This model simply adds a new contribution to 
the energy density on the r.h.s. of the FE 0, which would dominate at a late epoch in the 
evolution of the universe. In our article we consider similar modifications but, as we said, 
we constrain them in the BBN period. 

With all this in mind, we add a new term to Eq. ((7j) 



Po + A 



. y 



0.1 MeV 



(8) 



with 7 and A two arbitrary parameters, that we shall bound. Even if the type of modification 
that we introduce in (jHJ is arbitrary, some values of the parameter 7 have a particular 
physical significance. For example, the contribution of a term with 7 = would correspond 
to a cosmological constant. For 7 = 4, the new term would be of the form of extra relativistic 
particles, and for 7 = 3 it would correspond to extra massive particles. 

The new term in the FE modifies the expansion rate, which plays a crucial role in the 
dynamics during BBN. The universe, just before this period, consists of a hot dense plasma 
made of n, p, e^, v and 7. Neutrons and protons are in equilibrium through weak interac- 
tions. This equilibrium is lost when the universe is expanding so fast that particles cannot 
interact with each other. It happens when the expansion rate H equals the weak interaction 
rate Tw, a process that is known as the freeze out of the weak interaction. In general, any 
modification of the expansion rate changes the time of the freeze out. 

After freeze out, the temperature of the cooling universe will eventually be low enough 
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to permit the formation of the hghtest nuchdes. In order to calculate the abundances of 
these nuclides is important to know how many neutrons and protons are there when the 
nuclear reactions begin. The value n/p has its equilibrium value, which decreases with time 
{rrin > mp) , until the freeze out of weak interactions. Since then, neutrons are not destroyed 
or created by two body reactions, but can still be destroyed by neutron decay. Anyway, 
to know how many neutrons are there, one has to know at which temperature the weak 
reactions freeze out, because depending on the freeze out temperature, there will be more 
or less neutrons available to form the lightest elements. 

Before we do the numerical analysis, we can understand some qualitative features of the 
effect of the new term in (jSI). We can see that for positive values of A, its effect is to increase 
the expansion rate of the universe. Then, as we have said, the weak interaction will freeze 
out before the time of the freeze out in the standard case, and there will be more neutrons 
present in the universe at the moment the nuclides begin to form. So, the presence of a term 
with A > will increase the light element yields . 

The effect of a A < term would be the contrary, but in that case there is the problem 
that for certain values of A and 7, becomes negative, which is physically unacceptable. 
This is indeed what happens. If we look at (jHl), we can see that for large values of 7 (7 > 4) 
this occurs at high values of the temperature, i.e. when BBN begins (T ~ 1 — 10 MeV). On 
the other hand, for low or negative values of 7, becomes negative at low temperature, 
i.e. when BBN finishes (T ~ IkeV). Anyway, for all values of 7, there are problems when 
A < 0. We could be able to study particular cases of 7 where A < but, since we are 
interested in a systematic study, we will concentrate only on positive values of A. 

We have modified the standard numerical code jOI, performing the modification (jHJ and 
we have obtained the new predicted values for the several light element yields. An input 
that is needed for our calculation is the baryon-to-photon number ratio r]B = n^/n^. We 
use values consistent with the analysis of the cosmic microwave background anisotropies, 
where the parameter rjB plays a role [21, 

In Figs. [T(a)1and [I(b)j we can see the results for the deuterium and '^He abundances as 
a function of A, for several values of 7. The horizontal line corresponds to the observational 
upper limit |,8]. We are presenting the results for deuterium and "^He because these two 
elements constrain more A and 7 than the other elements. We can see that the abundances 
grow with A, in agreement with what we have reasoned before. 



In Fig. 121 the observationally allowed region for A and 7 is presented. Solid, dotted, and 
dashed lines correspond to the constraints coming from ^He, deuterium, and ^Li, respectively. 
The intersection of the constraints, giving the allowed region is shadowed. From this figure, 
we can see that the constraint coming from ^He is important for larger values of 7 and the 
one coming from deuterium is more relevant for smaller values of 7. We can understand this 
feature in the following way. When the nuclear reactions that produce the light elements 
begin, the first element formed from neutrons and protons is deuterium. Through ^He 
and tritium, all the deuterium converts into ^He, which is the most stable nuclide. It is 
when the reactions that convert deuterium into ^He freeze out that primordial deuterium 
begins to form. Therefore, primordial deuterium is formed after all the primordial ^He is 
formed. For large values of 7, the modification AT''' we have introduced is important at the 
beginning of BBN (high temperature), i.e., during the formation of ^He. For small values 
of 7, the modification is more important at the end of BBN (low temperature), i.e., during 
the formation of deuterium (when all the ^He is already formed). This feature can also be 
seen m Fig. IB In Fig. [T(I)| the maximum allowed value of A, given by the observational 
constraint from deuterium, is smaller for low values of 7 and grows as 7 becomes larger. 
The case of '^He, Fig. |l(b)[ is similar for the low values of 7, but for 7 > 3 the effect of the 
term AT'^ is more important, and the maximum allowed value of A is more constrained. 

From Fig. |21 we see that there are strong bounds on A, so we conclude that depar- 
tures from the FE are very much constrained. To weight the importance of the new term 
A (T/0.1 MeV)^ in (jH} we can compare it numerically with the standard contribution po in 
(0). For the maximum allowed values of A, and making the comparison at the standard 
BBN temperature T = 0.1 MeV, we find that the new term does not exceed 10^^ times the 
standard term. 

Carroll and Kaplinghat |.9j have done a similar study, motivated by the same concern 
than we have discussed in this section. They use the following generalization of the FE 



and bound the parameters a and Hi using BBN. The form Q is suitable to test radical 
departures from the FE as suggested by some work in brane- world models However, 
brane-world models developed in a series of papers that come closer to the standard FE 
are, in our opinion, more promising. Thus, we think that, in this sense, our parameteri- 




(9) 
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zation is more suitable. More support to our point of view and the form we use for the 
modified FE is the fact that, as it follows from our study, possibles departures from the 
FE have to be small due to the observational data on BBN. Even if both parameterizations 
are not mathematically completely equivalent, one may establish connections between the 
parameterization in Eq. Q and ours, in Eq. ©, when using suitable expansions in some 
limits. In this sense, our work complements the study in j^. 

Another related work has recently been done by Zahn and Zaldarriaga Q|. They show 
how future cosmic microwave background radiation experiments will probe the FE If 
we take their modification and test it in the radiation era using BBN (it would correspond 
to the particular case 7 = 4 in our Eg. ©) we reach bounds that are on the same order of 
magnitude that the values found in 10|. 



A Particular Case : Non-Universal Gravity 



In a recent paper [3] , Rafelski discusses a model where one would expect the fine structure 
constant a to vary in time. His work has been prompted by the recent measurements of 
a at cosmological time scales that may show this time variation effect He proposes 

that gravity couples much strongly (a factor ~ 10'^) to the second and third generation 
than it does to the first. In jsj it is claimed that the model is not excluded, since there 
is a lack of laboratory experimental constraints on how gravity couples to particles of the 
second and third generation (in the laboratory one measures the value of the gravitational 
strength through observations with matter of the first generation, namely electrons, protons 
and neutrons.) 

However, the claim in |2] is partially not true. While certainly the laboratory constraints 
on how muons, taus, and second and third generation quarks couple to gravity are practically 
absent, surprisingly enough there are strong constraints on the neutrino sector. Gasperini 



12| demonstrated that if gravity is not universally coupled to leptonic flavors, the gravi- 
tational field may contribute to neutrino oscillations. Laboratory experiments on neutrino 
osci 
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lations place very strong bounds on the strength of the violation of universal gravity 

i3- 



BBN also offers a test of this exotic possibility, since muonic and tauonic neutrinos 
participate actively in the expansion of the early universe. The FE in a model of non- 
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universal gravity reads 

H^ = ^[G (p^ + Pe + PuJ + G2 Pu, + Gs Pu.] (10) 

where G2 and G3 are the gravitational constants for the second and third generation and 
G is the usual Newton gravitational constant. As before, we have neglected the curvature 
term, which is not important in the radiation dominated epoch. 

We will now show that with some approximations we can use the general results that 
we have presented before in this section. First, we consider non-degenerate neutrinos (in 
the next section we shall drop this assumption). In this case, and if we ignore the small 
temperature difference between z/g, i/^, and we have p^^ = p^^ = p^^ = p^, with 

Pu = --T^ (11) 
^ 8 30 ^ ^ 



Then we can write ffTUI) as 

stt/^ r An^ 4- An„ 1 

(12) 



^2_8rfG 



AG2 + AGs 

Po + ^ Pu 



where po is given in (j3)) and we have defined AG2 and AGs as AGj = Gi — G. Since here 
we are interested in getting an approximate bound on AG2 and AG3, we consider that the 
temperature of the neutrinos and the photons are of the same order. Thus we can see that 
Eq. ()12j) is a particular case of the general modification (jH} with 7 = 4 and 

,.Y-^^^(,,U.Vf (13) 

Now, from Fig. |2l in the case 7 = 4, we can estimate the limit A < 10~^ MeV^, which 
constrains the deviation from universal gravity to values of the order 

< 1 (14) 

But before comparing this result with what is needed in ^j, we will see in the next section 
that if we allow for a neutrino degeneracy, all these constraints become weaker. 

III. TESTING THE FRIEDMANN EQUATION WITH DEGENERATE BBN 

Up to this point we have taken p^^ = 0, so we have not considered a neutrino/antineutrino 
asymmetry. In this section, we would like to study how the results change in the case of 
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degenerate neutrinos, i.e., ^ 0. To reduce the number of parameters, we put = fi^^ = 
fij,^ = fi^, namely that the asymmetry between neutrinos and antineutrinos is the same for 
the three famihes. Actually, it has been recently pointed out 14] that this condition is 
a consequence of oscillations among the different families of neutrinos. Because of the 
dependence of fi^ on the temperature, it is convenient to use the ratio = n^/Ty which is a 
constant quantity (again, we ignore a small temperature difference between and v^, Ur). 

The introduction of a neutrino degeneracy strongly affects the production of ^He. This 
can be easily understood because for this element, following a standard qualitative analysis 
one can get an approximate analytic expression for Yp, the primordial ^He mass fraction. 
It reads 

where {n/p)j^^^ is the ratio between neutrons and protons at the moment the ^He is formed. 
In the case of degenerate neutrinos this ratio is given by |3| 

(^/P) Nuc oc exp ( - -iu] (16) 

V Nuc J 

where we can see that degeneracy affects Yp exponentially. From Eqs. ()15|) and ()16|1 we can 
see that positive values of S,u diminish the predicted abundance Yp, an effect that goes in 
the opposite direction compared to the consequences of the modification (jH]), when A > 0. 
This can be seen in Fig. |3] where a value of = 0.06 increases the maximum allowed value 
of A by about one order of magnitude (such value of is consistent with recent analysis, 
see jlj]). As expected from our discussion above, the allowed region coming from the ^He 
restriction is significatively enlarged, while the restrictions coming from the other elements 
remain nearly constant. We conclude that for 7^ 0, larger values of A are allowed, so one 
has much more freedom to modificate Eq. (0) than in the case = 0. 



Constraints on Non-Universal Gravity from Degenerate BBN 

In the same way as we did in the precedent section, we could treat the model of non- 
universal gravity as a particular case of (jH)), and from Fig. El we could limit the parameter 
AG2+AG3 ^Q], j^Yie value = 0.06. But here, taking profit that in this case we have fixed one 
parameter (7 = 4), we can make a more general analysis of how affects the restrictions 
on A. 
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We notice that Eq. (fT^ is still valid in the degenerate case (recall we set fi^^ = jiy^ = 
fi^^ = fi^) but now, to get p^, we have to use Eq. ® with the temperature T equal to 
the neutrino temperature, as it should be. We make a numerical analysis with the free 
parameters ^'^a+AGa -^^ -p[g, |^we show the allowed region for the parameter space 

^ AG2+AG3 ^^^^ for the value r]B = 4: x 10^^°. We can see that constraints coming from the 
observational limits of the different elements bound these parameters into a closed region. 
Notice that if one only considered the restrictions coming from one element, the allowed 
region would be no longer closed. The reason for this is that the effects of '^'^a+^'^s and 

on each element abundance go in opposite directions and may compensate each other. 
As we have noted before in the case of ^He, positive values of diminish the predicted 
abundance Y.p. On the other hand, positives values of ^"^z+ags increase the expansion rate 
of the universe, H, which leads to an increase of Yp, compensating the change induced by 

In Fig. El we show the allowed region for different values oi rjB- We note that for any 
r]B-i there are values of ^'^2 +^'^3 g^^^ ^j-^g^^^ gjyg theoretical predictions compatible with 
the primordial element observations. Rather than give accurate bounds on ( ^'^^g^'^^ ,^^), 
here we are interested in pointing out that the deviation from universal gravity is roughly 
limited to values I AG2+AG3 1 < 20. At the light of this result, we conclude that the values of 
G2 and 6*3 that are postulated in [3|, i.e., about a factor of 1,000 bigger than the value of 
the Newton gravitational constant (for the first family), are excluded. 

Our bounds on non-universal gravity are much weaker than the bounds coming from 

n 

oscillation experiments, which are at the level ~ 10 jl^. However, the bounds presented 
in this work may have some interest due to the following reasons. First of all, our limits 
and the limits from oscillations are obtained at very different times; BBN period {z ~ 10^) 
versus today {z = 0). Also, there might be instances where oscillation experiments do 
not bound non-universal gravity. Two such cases are discussed by Gasperini Q]. First, 
if there is no mixing between neutrino weak eigenstates and the total energy eigenstates 
there is no gravitational induced oscillations, and thus no bounds. Second, there might be 
an anomalous gravitational coupling of the neutrino with respect to the charged fields but 
with the three u coupling equally, namely, I'ej^fi, and u^- all couple with strength G' 7^ G. 
Then again there are no gravitationally induced oscillations. However, our limits still apply 
in these two cases. 
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FIG. 1: Yields of primordial deuterium, normalized to hydrogen, (a), and mass fraction of primor- 
dial helium (b), in a universe with evolution in the BBN epoch given by Eq. ((HJ, as a function of A 
and for 7 = —3, 0, 3, and 5. We fix r/^ = 4 x 10^^''. The horizontal line of each figure corresponds 
to observational data. We only show the upper limit, since the lower limit is below the displayed 
scale. Thus, the region above the horizontal line is excluded. 
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FIG. 2: Constraints on the (7, A) plane coming from the requirement that the predicted abundances 
are consistent with the observed values, in the case fii, = 0. We show the constraints from "^He as 
a solid line, from D as a short-dashed line, and from ^Li as a long-dashed line. The intersection 
gives the allowed region, that we shadow. We have put 7?b = 4 x 10^^*^. 




FIG. 3: Same that Fig. 21 but with S^i, = Hu/T = 0.06. The constraints are much looser as can be 
seen by noticing the scale change in the A axis. 
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FIG. 4: Region in the ^ plane where the modified Friedmann equation leads to 

abundances in agreement with observation. We display the limits coming from ^He (solid line), D 
(short-dashed line), and ^Li (long-dashed line). The allowed region when one takes into account 
all constraints is shadowed. In this figure we set ry^ = 4 x 10"^*^. 
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FIG. 5: Allowed region in the ( ^'^^g^'^^ ; ^v) plane. We vary r\B'- we show as a solid line the case 
r]B = 4: X 10^^'^ (same that in Fig. 0]), and as a short-dashed line and as a long-dashed line the 
cases r]B = 3 X 10^^^ and rjB = ^ x 10^^*^, respectively. 
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